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ABSTRACT 


Further  analysis  of  the  anomalous  events  in  the  eastern  Himalaya  region 
tfith  an  expanded  set  of  stations  shows  that  the  Ms  values  of  these  events  are 
indeed  lew  compared  to  ny  Other  characteristics  of  these  events  indicate, 
however,  that  all  these  events  are  earthquakes  rather  than  explosions. 
Dilatational  first  motion  for  P waves,  long  period  S to  Rayleigh,  short 
period  S to  P wave  amplitude  ratios  are  characteristic  of  earthquakes. 
Readings  of  the  PP  phase  where  available  as  well  as  epicenter  calculations 
indicate  that  the  anomalous  event  hypocenter  are  shallow,  i.e.,  less  than 
80  km.  Data  quality  did  not  permit  the  refinement  of  epicenter  locations. 

A search  for  additional  anomalous  events  with  respect  to  M -nr  in  a 
wider  area  in  the  Himalayan-Tibet  region  indicated  that  there  are  no  other 

contiguous  areas  containing  anomalous  events  in  this  region  beyond  those 
already  found. 
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INTRODUCTION 


This  report  is  a supplement  to  a previous  report  on  low  M events  in  the 
Himalayan  region  (Der,  1973),  hereafter  referred  to  as  Report  //296.  It 
includes  additional  data  for  a larger  network  of  stations,  an  evaluation  of 
discriminants  other  than  Mg  vs.  m^,  and  a search  for  low  Mg  events. 

We  first  reanalyze  the  anomalous  events  from  296  with  respect  to  depth 
as  determined  from  P and  S-P  times,  and  from  P-pP  times.  Discrimination  by 
means  of  short-period  S to  P amplitude  ratios  is  discussed,  and  the  possible 
use  of  LQ/LR  ratios  is  explored.  Then  the  effects  on  Mg— m,  values  of  more 
complete  P wave  data  is  explored  together  with  the  expected  effects  on  M of 
an  abnormally  thick  crust.  Finally  we  report  on  the  use  of  station  KBL  to 
search  for  anomalous  events  in  a wider  area  than  searched  in  SDL  296. 


DATA 


The  events  selected  from  Report  if 296  for  this  study  are  shown  in  Table  I. 
TVo  criteria  were  used  in  the  selection:  the  low  Ms  value  compared  to  n^,  and 

the  proximity  of  the  epicenter  to  30°N  and  95°E  regardless  of  relative  Mg-mb 
values,  although  most  of  the  30°N,  95°E  events  also  qualify  as  'Mow  M " S 

g 

earthquakes . 

In  addition  to  the  station  network  used  in  the  previous  study  (Der,  1973) 
we  acquired  data  from  the  WWSSN  stations  AAE,  BUL,  CTA,  COL,  SEO,  TRI  and  NUR. 
The  map  on  Figure  1 shows  the  location  of  all  the  stations  used.  The  basic 
data  consists  of  short  and  long  period  seismograms  on  35  mm  film  chips.  The 
quality  of  recordings  is  often  poor,  and  the  quality  of  stations  is  often 
very  unequal.  Some  operate  at  high  gains  and  are  very  useful  (KBL  for 
example)  while  some  (LAH  for  example)  are  almost  useless.  Most  of  the 

recordings  are  not  suitable  for  digital  conversion  because  of  the  poor  data 
quality. 
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TABLE  I 

LIST  OF  EVENTS  STUDIED 


Yr 

Mo 

Day 

Hr 

Min 

Sec 

Coordinates 
N.  Lat  E.  Long 

SDL 

M 

— s 

296 

2b 

New 

Avg. 

"b 

Depth 

1 

64 

10 

06 

2 

54 

32.7 

30.3 

94.6 

3.61 

- 

4.54 

N* 

2 

66 

07 

05 

10 

1 

22.0 

27.5 

92.4 

3.90 

5.32 

5.32 

77 

3 

66 

09 

11 

15 

55 

20.0 

27.0 

95.8 

3.77 

5.24 

4.69 

37 

4 

66 

09 

26 

6 

3 

48.0 

27.6 

92.7 

3.59 

4.49 

5.39 

N 

5 

67 

07 

07 

22 

56 

30.8 

27.8 

92.2 

3.71 

4.76 

4.78 

f. 

6 

68 

06 

28 

20 

34 

55.3 

30.1 

95.1 

3.54 

4.77 

4.82 

44 

7 

68 

06 

30 

5 

4 

10.0 

30.2 

94.8 

3.40 

4.61 

4.64 

42 

8 

68 

07 

01 

3 

11 

10.0 

30.3 

94.5 

3.00 

4.14 

4.40 

28 

9 

68 

07 

04 

6 

45 

58.0 

30.3 

94.9 

3.46 

4.87 

4.82 

N 

10 

68 

07 

13 

6 

5 

54.2 

30.3 

94.6 

3.46 

4.73 

4.75 

N 

11 

68 

07 

14 

18 

12 

41.0 

30.3 

94.8 

3.54 

4.72 

4.89 

22 

12 

68 

07 

16 

22 

23 

7.0 

30.3 

94.8 

3.45 

4.57 

4.70 

40 

13 

68 

07 

19 

18 

48 

59.0 

30.2 

94.9 

3.56 

- 

4.80 

N 

14 

68 

07 

23 

20 

51 

47.9 

30.3 

94.9 

3.43 

4.57 

4.69 

30 

15 

68 

07 

25 

3 

34 

13.0 

30.2 

94.8 

3.27 

- 

4.90 

N 

16 

68 

07 

26 

12 

44 

3.0 

29.4 

95.0 

3.45 

4.65 

4.73 

N 

17 

68 

08 

23 

12 

1 

16.5 

30.3 

94.9 

3.46 

4.62 

4.75 

N 

18 

68 

08 

25 

17 

55 

5.3 

30.4 

94.8 

3.28 

4.39 

4.58 

19 

19 

68 

08 

29 

19 

51 

24.6 

30.2 

95.1 

3.48 

4.55 

4.68 

N 

20 

68 

09 

01 

5 

59 

26.6 

30.3 

94.8 

3.59 

4.56 

4.60 

20 

21 

68 

09 

03 

17 

45 

54.1 

30.2 

94.8 

3.43 

4.52 

4.64 

53 

22 

69 

08 

15 

7 

15 

37.0 

30.2 

95.0 

3.57 

4.87 

4.90 

N 

N 


normal  depth 


33  km 


RESULTS  OF  DATA  ANALYSIS 


SHORT- PERIOD  WAVE  TRAVEL  TIMES 

Table  II  shows  the  arrival  timer,  of  short  period  P and  S waves  read  from 
the  iilm  chips.  Arrival  times  with  estimated  uncertainties  larger  than  3 
sec  were  omitted  from  the  table. 

P WAVE  FIRST  MOTIONS 

No  clear  long-period  P waves  were  found  for  the  clustered  events,  but 
short-period  P wave  first  motions  have  been  observed  for  many  of  the  events. 
The  readings  and  their  quality  classifications  are  given  in  Table  HI. 

Single  events  do  not  have  enough  readings  to  define  a fault  plane  solution. 
Nevertheless  clear  dilatations  have  been  observed  (Figure  2a).  We  considered 
the  possibility  that  the  events  clustered  at  35°N,  95°E  could  possess  a 
common  source  mechanism  as  many  earthquake  swarms  and  aftershock  sequences 
do.  The  combined  first  motions,  however,  for  the  event  cluster  do  not  show 
a pattern  consistent  with  any  common  fault  plane  solution.  The  presence  of 
dilatations  indicates,  however,  that  the  events  in  question  are  earthquakes 
and  not  explosions. 

SHORT-PERIOD  WAVE  AMPLITUDES 

Table  IV  shows  the  amplitudes  and  dominant  periods  of  short  period  P 
and  S waves  at  stations  where  both  types  of  waves  could  be  observed. 

Von  Seggern  (1972)  showed  that  the  ratio  of  the  amplitudes  of  short-period 
P and  S waves  is  a discriminant  between  earthquakes  and  explosions.  The 
ratios  in  Table  IV  are  characteristic  of  earthquakes  as  seen  by  comparison 
with  Figure  2b  taken  from  the  report  by  von  Seggern. 

pF  PHASES 

The  pP-P  times  read  at  teleseismic  stations  are  shown  in  Table  V.  The 
pP  phases  are  uncertain  and  the  times  do  not  have  the  proper  moveout  for 
pP  phases,  this  indicates  that  some  of  the  phases  were  incorrectly  identified. 
The  depth  range  consistent  with  the  pP-P  times  are  also  given  in  the  table. 

All  the  depths  are  shallow,  and  are  not  inconsistent  with  National  Ocean 
Survey  (NOS)  depths  in  Table  I. 
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TABLE  II 

Travel  Times  of  P and  S Waves 


Event 


No. 

Date/Time 

Station 

P-Tirae 

S-Time 

1 

64/10/06 

CHG 

02:57:24.7 

t=02:54:33.0 

COL 

03:06:03.5 

NDI 

02:58:08.0 

SHI 

03:01:30.6 

SHL 

02:55:52.0 

TRI 

03:04:55.3 

2 

66/07/05 

NDI 

10:04:28.8 

10:06:55.0 

t=10:01:22.0 

QUE 

10:06:19.0 

10:10:33.6 

SHL 

10:01:55.2 

3 

66/07/05 

NDI 

16:07:30.0 

t=15  :55 :20 

CHG 

15:57:29.5 

15:57:29.5 

COL 

16:07:05.0 

NDI 

15:59:04.5 

16:01:55.2 

NUR 

16:05:05.0 

SEO 

16:01:09.3 

SHI 

16:02:35.2 

SHL 

15:56:18.0 

4 

66/09/26 

CHG 

06:06:17.3 

COL 

06:15:39.0 

NDI 

06:07:00.7 

POO 

06:08:16.0 

QUE 

06:08:50.0 

SHL 

06:04:23.3 

5 

67/07/07 

BUL 

23:08:28.2 

t=22 :56 : 30. 8 

CHG 

22:49:06.5 

COL 

23:08:22.0 

NDI 

22:59:30.6 

23:01:51.0 

QUE 

23:01:25.7 

23:05:41.0 

SHI 

23:03:19.0 

TRI 

23:06:54.5 

6 

68/06/28 

BUL 

20:47:08.0 

t=20:34:55 

CHG 

20:37:47.0 

COL 

20:46:25.2 

CTA 

20:46:03.0 

KBL 

20:39:49.2 

20:43:53.0 

QUE 

20:40:10.5 

SHL 

20:36:14.9 

TRI 

20:45:20.0 

Component 


SPE 
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TABLE  II  (Continued) 

Travel  Times  of  P and  S Waves 


Event 


No. 

7 


10 


11 


12 


Date/ Time 

Station 

P-Time 

S-Titne 

68/06/30 

BUL 

05:16:22.8 

t=05:04:10.0 

CHG 

05:07:02.6 

COL 

05:15:39.5 

KBL 

05:09:04.1 

05:13:08.1 

SHL 

05:05:30.0 

05:06:21.5 

68/07/01 

BUL 

03:23:20.0 

t=03 : 11 :10.0 

CHG 

03:14:02.7 

COL 

03:22:37.0 

SHL 

03:12:28.9 

03:13:26.0 

68/07/04 

BUL 

06:58:11.0 

t=06 :45 :58.0 

COL 

06:57:29.0 

KBL 

06:50:52.3 

06:54:56.1 

NUR 

06:55:21.0 

SHI 

06:53:02.9 

68/07/13 

BUL 

06:18:05.0 

CHG 

06:08:45.8 

COL 

06:17:22.5 

CTA 

06:17:07.2 

KBL 

06:10:47.6 

06:14:51.0 

NDI 

06:09:26.2 

NUR 

06:15:19.5 

QUE 

06:11:04.5 

SHI 

06:12:56.5 

TRI 

06:16:17.0 

68/07/14 

BUL 

18:24:56.2 

t«18:21:41 

CHG 

18:15:35.3 

COL 

18:24:11.0 

KBL 

18:17:36.4 

18:21:40.1 

NDI 

18:16:15.5 

NUR 

18:22:07.5 

POO 

18:17:40.0 

SHI 

18:19:45.5 

SHL 

18:14:02.8 

TRI 

18:23:08.5 

68/07/16 

BUL 

22:34:19.5 

t«22:23:07.0 

CHG 

22:25:56.2 

COL 

22:34:36.0 

CTA 

22:34:17.4 

KBL 

22:28:00.9 

22:32:04.0 

"V 


Component 


SPN 


SPN 


SPN 


SPE 


SPE 


SPN 


'X 
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TABLE  II  (Continued) 


Event 

No. Tate/ Time 


13  68/07/19 

t= 1 8:48: 59 


14  68/07/23 


15  68/07/25 

t=03:34:13 


16  68/07/26 

t=12 :44:03 


Travel  Times  of  P and  S Waves 


Station 

P-Time 

S-Time 

Component 

NDI 

22:26:39.3 

NUR 

22:32:32.2 

SHI 

22:30:09.9 

SHL 

22:24:26.9 

AAE 

18:58:41.5 

B'JL 

19:01:11.9 

CHG 

18:51:51.8 

COL 

19:00:30.7 

KBL 

18:53:54.5 

18:57:57.0 

SPN 

NDI 

18:52:32.4 

NUR 

18:58:25.9 

QUE 

18:54:12.8 

18:58:33.0 

SPN 

SHI 

18:56:03.0 

SHL 

18:50:19.5 

TRI 

18:59:24.9 

BUL 

21:04:02.0 

CHG 

20-54:41.0 

COL 

2.1  :03:18. 5 

KBL 

20:56:43.4 

21:00:46.4 

SPN 

NDI 

20:55:21.5 

NUR 

21:01:13.0 

SHL 

20:53:08.5 

SHI 

20:58:50.8 

BUL 

03:46:25.5 

CHG 

03:37:06.2 

COL 

03:45:43.0 

NDI 

03:37:46.8 

NUR 

03:43:37.0 

SHI 

03:41:15.5 

AAE 

12:53:44.0 

BUL 

12:56:14.5 

CHG 

12:46:59.5 

COL 

12:55:37.9 

CTA 

12:55:15.0 

KBL 

12:49:02.0 

12:53:05.0 

SPE 

NUR 

12:53:33.5 

QUE 

12:49:23.3 

SHI 

12:51:11.0 

SHL 

12:45:12.7 

12:46:12.5 

SPS 

TRI 

12:54:38.0 
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TABLE  II  (Continued) 

Travel  Times  of  P and  S Waves 


Event 


68/08/23 
t=12 :01 :16 


B'JL 

CHG 

COL 

CTA 

KBL 

NDI 

NUR 

QUE 

SHL 


12:13 

12:04 

12:12 

12:12 

12:06 

12:04 

12:10 

12:06 

12:02 


34.5 

07.7 

45.0 

31.8 

10.7 

50.0 

42.5 

35.7 

37.0 


12:10:14.1 


12:03:34.2 


SPN 


SPE 


18  68/08/25 

t=17:55:05 


19  68/08/29 

t-19 :51 : 25 


20  68/09/01 

t=05:59 :27 


21  68/09/03 

t«*17:45:54 


BUL 

18:07:20.0 

CHG 

17:57:59.4 

COL 

18:06:38.0 

CTA 

18:06:19.3 

KBL 

18:00:01.4 

NDI 

17:58:41.1 

NUR 

18:04:33.0 

QUE 

18:00:20.2 

SHL 

17:56:27.5 

BUL 

20:03:37.5 

CHG 

19:54:17.3 

COL 

20:02:55.0 

CTA 

20:02:35.0 

KBL 

19:56:19.0 

NDI 

19:54:59.4 

NUR 

20:00:49.7 

SHI 

19:58:28.8 

SHL 

19:52:45.0 

TRI 

20:01:50.0 

BUL 

06:11:41.5 

CHG 

06:02:19.5 

COL 

06:10:59.0 

CTA 

06:10:37.9 

KBL 

06:06:06.0 

SHI 

06:06:32.0 

SHL 

06:00:48.4 

CHG 

17:48:45.0 

COL 

17:57:22.0 

CTA 

17:57:03.0 

KBL 

17:50:47.2 

18:04:04.5  SPE 


17:57:36.7  SPE 


20:00:22.8  SPN 


06:08:26.0  SPN 


17:54:50.5  SPE 
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TABLE  II  (Continued) 

Travel  Times  of  P and  S Waves 


Event 
Date/ Time 


69/08/15 
t=07 : 1 5 : 37 


Station 

P-Time 

S-Time 

NDI 

17:49:26.0 

NUR 

17:55:16.5 

POO 

17:50:49.3 

QUE 

17:51:08.7 

SHL 

17:47:12.1 

17:48:08.0 

CMC 

07:18:42.8 

COL 

07:27:07.9 

KBL 

07:20:32.7 

07:24:37.0 

NIL 

07:19:51.6 

07:23:27.0 

QUE 

07:20:52.5 

SHI 

07:22:42.0 

SHL 

07:16:57.4 

07:17:51.7 

Component 


! 


TIME  MARK 


SHL 


16  FEB  68 


START  / 

If  j /period  too\ 


TIME  MARK 

I 


START 


TIME  MARK 


START 


START 


SHl  ^ 1, 

V“'-v  ^ (unreadable) 

v6  JUN  68  \ 

START  TIME  MARKS 
CHG  | || 


KBL 


TIME  MARK 

n 


v ^ 

21  AUC  68 


SHL 


14  JUL  68 


START 
TIME  MARK 


,4 


KBL 


START 

V*^A/VA^/\/^ 

23  JUL  68  y 


START 


SHL 
k/-V 
23  AUG  88 


SHL 


01  SEP 


START 

« ~^jl 


TIME  MARK 

n 


START 

03  SEP  86  * f 


TIME  MARK 
.V^A 


'V' 


START 

11  IIP  68  lfni| 


TIME  MARK 


TIME  MARK 


START 


SHL 

*y  "*■  • ~»nv\j< 


TIME  MARK 

l 


18  SEP  68 

Figure  2a.  Dilatational  first  '.otion  observed  for  various  even's. 

-18- 


I 


'V. 


*r"-T' 


tabu:  iv 


Short  Period  S/P  Amplitude  Ratios 


Event 

Stat ion 

R(S/P  Amplitude 
Ratio) 

Log  R 

7/0  S/66 

QUE 

1.4 

0.16 

NDI 

28.4 

1.45 

7/07/67 

QUE 

5.5 

0.74 

NDI 

29.7 

1.47 

6/28/68 

KBL 

1 .3 

0.10 

6/30/68 

CHG 

2.6 

0.41 

SHL 

16.0 

1 . 20 

KBL 

0.4 

-0.45 

7/01/68 

SHL 

77.1 

1.89 

7/26/68 

KBL 

1.1 

0.05 

SHL 

65.4 

1.82 

8/23/68 

SHL 

31  . 1 

1 .49 

8/25/68 

SHL 

1.1 

0.06 

8/29/68 

SHL 

0.8 

-0.11 

9/03/68 

SHL 

15.0 

1.18 

8/15/69 

SHL 

7.8 

0.89 

Tt  ELE  V 
pP  - P Times 


SEQ 

SHI 

NUR 

AAE 

66/09/11 

- 

_ 

_ 

67/07/07 

- 

5.5 

_ 

68/06/28 

- 

- 

68/07/01 

- 

NR 

NR 

_ 

68/07/13 

- 

4.6 

+ 

68/07/14 

7 

5.0 

5.0 

68/07/16 

- 

10.1 

+ 

_ 

68/07/23 

- 

7.4 

7 

68/07/25 

- 

9.5 

+ 

68/07/26 

- 

10.5? 

+ 

+ 

68/08/25 

- 

NR 

5.0 

68/08/29 

- 

7.9 

+ 

68/09/01 

7 

+ 

— 

_ 

69/08/15 

- 

+ 

? 

- 

+ P Observed,  no 

PP 

- No  P 
NR  No  Record 
• Doubtful  Readings 


CTA 

COL 

BUL 

TRI 

Avg 

DEPTHfL'm'* 

- 

6.0 

- 

20 

- 

+ 

+ 

+ 

17 

5.5 

+ 

+ 

+ 

17 

- 

5.0 

+ 

- 

16 

+ 

7.5 

+ 

2.5 

10-24 

- 

6.5 

+ 

11.5 

17-40 

+ 

9.5 

+ 

- 

-35 

- 

+ 

5.9 

- 

15-27 

- 

+ 

+ 

- 

35-40 

5.0 

+ 

+ 

4.5 

14-15 

4.2 

+ 

+ 

- 

13-18 

+ 

+ 

5.0 

+ 

15-27 

+ 

+ 

4.5 

7 

12 

- 

10.1 

- 

- 

33 
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DESCRIPTION  OF  LONG-i  ERIOD  SEISMOGRAMS 
Events  Clustered  at  30°N  and  95CE 


Figures  I through  9 show  long-period  three-component  seismograms  at 
Shillong  (SHL)  for  seven  events.  SHL  is  the  only  station  which  recorded  the 
surtace  waves  well  from  the  event  cluster  at  30°N,  95°E.  Weak  surface  waves 
were  observed  at  other  stations,  especially  at  Chiang  Mai  (CHG) , but  these 
are  hardly  worth  reproducing  due  to  the  low  signal-to-noise  ratio.  The  recordings 
at  Shillong  all  show  a Rayleigh  pulse  on  the  vertical  component  which  regains 
reasonably  similar  from  event  to  event.  The  group  velocity  of  the  Rayleigh 
pulse  is  about  3km/sec  and  the  period  is  about  10  sec.  The  horizontal  com- 
ponents are  also  very  similar  from  event  to  event  to  a very  fine  level  of 
detail.  Short-period  waves  corresponding  to  Lg  arrive  before  the  Rayleigh 
pulse  and  there  are  indications  of  longer  period  waves,  probably  fundamental 
Love,  preceding  the  Rayleigh  pulse. 


The  Rayleigh  pulse  iv  also  seen  on  the  horizontal  components. 


me  i fiLi 

Of  the  horizontal  to  vertical  amplitude  is  close  to  .7  ruling  out  the  possi- 
bility of  the  pulse  consisting  of  higher  modes  since  for  higher  modes  this 
ratio  should  be  much  lower  for  any  reasonable  mantle  structure.  The  wave- 
forms are  quite  similar,  suggesting  similar  source  mechanisms  for  these 
events. 


Discussion  of  Long-Period  Seismograms  of  Events  not  Located  at  30°N.  95°E 

The  most  prominent  feature  of  long-period  seismograms  from  the  10/06/64 
event  is  the  Rayleigh  wave  (Figure  10).  Higher  mode  surface  waves  may  be  also 
present  in  the  higher  frequency  waves  preceding  the  Rayleigh  pulse.  The 
quality  of  seismograms  is  poor. 

Figure  11  shows  seismograms  of  the  event  7/05/66  recorded  at  Shillong 
at  a distance  of  220  km.  P and  S phases  are  visible.  There  is  an  indication 
of  Love  waves.  The  depth  of  this  event  was  given  by  the  NOS  as  77  km. 

Figure  12  shows  seismograms  for  the  event  9/11/66.  At  Shillong  long- 
period  P and  S phases  are  faintly  visible.  Clear  Rayleigh  and  Love  waves 
are  present.  There  are  indications  of  higher  Rayleigh  modes  in  the  form  of 

short-period  waves  preceding  the  fundamental  Rayleigh  mode  on  the  vertical 
component. 


-22- 


At  Chiang  Mai  no  P or  S is  visible.  Fundamental  Rayleigh  and  Love  waves 
are  present,  as  well  as  indications  of  higher  mode  Rayleigh  waves. 

The  recordings  of  the  earthquake  9/26/66  at  Shillong  (Figure  13) 
again  show  long  period  P,  S,  and  Rayleigh  waves. 

The  recording  of  the  event  7/07/67  shows  S and  Rayleigh  waves  (Figure  14). 

Long-Period  S to  Long-Period  Rayleigh  Amplitude  Ratios 

Von  Seggern  (1972)  and  Blandford  and  Clark  (1974)  showed  that  the  long- 
period  S vs.  Rayleigh  wave  amplitude  ratio  can  be  a good  discriminant  between 
earthquakes  and  explosions.  Explosions  generate  long-period  shear  waves  much 
less  efficiently  than  earthquakes  if  the  events  are  scaled  to  the  same  maxi- 
mum Rayleigh  amplitude.  The  data  presented  by  von  Seggern  show  that  the  mean 
of  the  logarithm  of  S/Rayleigh  amplitude  ratios  is  around  -.3,  corresponding 
to  a ratio  of  about  .5.  Only  the  largest  explosions  generate  visible  long- 
period  S waves.  Since  most  of  the  events  investigated  in  this  report  show 
at  least  indications  of  long-period  S waves,  this  fact  by  itself  place  them 
in  the  earthquake  population.  Table  VI  shows  the  observed  S to  Rayleigh 
amplitude  ratios  at  SHL,  The  values  were  obtained  by  measuring  the 
periods  and  peak-to-peak  amplitudes  of  the  waves  and  correcting  the  latter 
for  the  instrument  response.  Rayleigh  wave  amplitudes  were  measured  on  the 
vertical  component,  S wave  amplitudes  on  the  horizontal  components  where  the 
amplitude  was  maximum. 

The  value  for  the  clustered  events  is  based  on  some  faint  indications  of 
long  period  S.  The  measurements  are  questionable  and  if  real  they  represent  a 
maximum  value.  The  ratios  fall  within  the  earthquake  region  of  von  Seggern, 
although  such  ratios  are  rather  crude  discriminants  in  the  absence  of  dis- 
tance corrections. 

Location  and  Depth  Studies 

Table  VII  shows  the  results  of  location  and  depth  determination  using  the 
travel  times  in  Table  II  and  with  no  constraints  on  either  depth  or  location. 
Most  of  the  depths  are  moderate  except  events  1 and  9;  some  are  negative. 


I 

i 

1 


I 

I 


TABLE  VI 

Long  Period/Ray leigh  Amplitude  Ratios 


Event 

R (S/Rayleigh 
Amplitude  Ratio) 

Log  : 

10/06/64 

.75  ? 

-.12 

07/05/66 

.25 

-.61 

07/07/66 

.62 

-.21 

09/11/66 

.60 

-.23 

09/26/66 

.78 

-.11 

Cluster  at 
30  °N,  95°E 

-.84  ?? 

-.08 

I 
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TABLE  VII 

Relocation  of  Events  Without  Constraints 


Mo 

Date 
Day  Yr 

Origin  Time 

Depth 

SDV(R) 

Shift 

kTM  A7 

10 

06 

64 

02:54:45.9 

175 

0.8 

92 

201 

09 

11 

66 

15:55:12.1 

- 34 

0.6 

59 

348 

09 

26 

66 

NO  CONVERGENCE 

07 

07 

67 

22:56:35.1 

49 

1.2 

16 

207 

06 

28 

68 

20:35:00.2 

74 

1.4 

17 

261 

06 

30 

68 

05:04:13.7 

52 

0.7 

7 

308 

07 

01 

68 

03:11:10.6 

46 

0.0 

21 

350 

07 

04 

68 

06:46:22.5 

227 

0.8 

78 

306 

07 

13 

68 

06:05:57.0 

50 

1.2 

5 

169 

07 

14 

68 

18:12:40.1 

1 

1.0 

24 

357 

07 

16 

68 

22:23:06.4 

23 

1.1 

6 

56 

07 

19 

68 

18:49:00.2 

16 

1.3 

23 

307 

07 

23 

68 

20:51:50.3 

31 

0.9 

19 

293 

07 

25 

68 

03:34:16.5 

49 

0.2 

21 

302 

07 

26 

68 

12:44:12.6 

67 

4.1 

48 

293 

08 

23 

68 

12:00:53.1 

146 

1.7 

56 

26 

08 

25 

68 

17:55:11.5 

54 

0.6 

24 

209 

08 

29 

68 

19:51:29.3 

62 

0.8 

33 

269 

09 

01 

68 

05:59:33.5 

70 

0.9 

25 

170 

09 

03 

68 

17:45:53.9 

32 

1.6 

18 

18 

08 

15 

69 

NO  CONVERGENCE 

As  the  next  step  we  tried  to  derive  regional  Lravel  time  anomalies  by 
constraining  the  depths  of  events  13,  14,  and  16  (which  were  recorded  at  many 
stations)  to  the  NOS  depths  (which  incidentally  are  within  the  range  allowed 
by  pP  times)  performing  a location  run,  and  considering  the  travel  time 
residuals  as  travel  time  anomalies  for  the  region  surrounding  these  events. 
This  is  the  SRST  approach  of  Veith  and  Clawson  (1973)  The  residuals  for  the 
three  events  showed  some  resemblance,  so  we  took  the  averages  for  the  three 
events  as  travel  time  anomalies.  These  are  shown  in  Table  VIII.  Applying 
the  residuals  as  travel  time  corrections  to  all  events,  we  relocated  the 
events.  The  results  of  the  computations  are  shown  in  Table  IX.  There  are 
no  significant  chai  ges  in  the  depths. 

An  attempt  was  also  made  to  refine  the  location  by  means  of  an  additional 
restraining  of  the  origin  time  by  use  of  S-P  arrival  time  differentials. 

S-P  times  together  with  the  standard  JB  tables  determined  a set  of  distance 
and  corresponding  P travel  times.  By  subtracting  the  travel  time  from  the 
arrival  time,  a origin  time  estimate  is  obtained.  However,  the  resulting 
origin  times  deduced  from  different  stations  generally  did  not  agree  within 
an  acceptable  limit.  The  reason  was  probably  in  the  difficulty  of  determining 
a good  S start  time  for  the  distant  stations.  Since  most  events  had  readings 
fiom  SHL  and  KBL,  at  distances  of  about  5°  and  20°  respectively,  it  was 
decided  to  use  only  these  stations  for  estimating  origin  times. 

A plot  of  the  observed  values  of  S-P  vs.  P time  for  those  events  with 
readings  at  both  SHL  and  KBL  (Figure  15) , indicates  that  the  average  slope 
is  about  .83  as  opposed  to  a slope  of  .78  from  the  JB  tables  within  the 
same  distance  range.  Therefore,  a constant  fraction,  .83/. 78,  of  the  JB 
S-P  versus  0-P  line  was  used  to  determine  th"  origin  time  from  readings  at 
SHL  and  KBL  for  all  events. 

Based  upon  the  new  origin  times,  location  runs  were  again  made  and  the 
results  are  shown  in  Tables  X and  XI.  As  can  be  seen  from  the  tables,  the 
depths  recalculated  do  not  change  appreciably.  Some  move  deeper,  some 
shallower  in  comparison  to  the  previous  runs,  but  with  the  exception  of  one 
event,  they  are  all  shallower  than  100  kms. 
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TABLE  VIII 

Travel  Time  Anomalies  (Residuals) 


SHL 

0.4 

CHG 

0.5 

NDI 

-0.9 

KBL 

1.3 

QUE 

-0,8 

SHI 

0.0 

NUR 

1.3 

TRI 

0.0 

COL 

-1.3 

BUL 

-0.9 

POO 

0.8 

SEO 

0.0 

0.0 

-0.4 


AAE 

CTA 


TABLE  IX 


Relocation  of  Events  with  Travel  Time 
Corrections  using  Residual 


Mo 

Date 
Day  Yr 

Origin  Time 

Depth 

SDV(R) 

Shift 
KM  A 7. 

10 

06  64 

02:54:43.5 

150 

0.8 

65 

203 

09 

11  66 

15:55:13.7 

- 23 

1.1 

59 

346 

09 

26  66 

NO  CONVERGENCE 

07 

07  67 

22:56:34.8 

50 

0.7 

12 

225 

06 

28  68 

20:35:01.3 

89 

1.3 

15 

247 

06 

30  68 

05:04:14.5 

67 

1.3 

3 

58 

07 

01  68 

03:11:10.3 

52 

0.0 

29 

348 

07 

04  68 

06:46:15.1 

176 

0.4 

23 

290 

07 

13  68 

06:05:58.0 

61 

2.0 

12 

187 

07 

14  68 

18:12:41.6 

11 

1.8 

24 

3 

07 

16  68 

22:23:08.3 

39 

1.8 

5 

93 

07 

19  68 

18:49:01.9 

29 

1.9 

30 

299 

07 

23  68 

20:51:51.5 

42 

1.8 

30 

95 

07 

25  68 

03:34:17.0 

56 

1.1 

23 

293 

07 

26  68 

12:44:14.0 

84 

4.2 

44 

284 

08 

23  68 

12:00:56.5 

-122 

2.0 

51 

26 

08 

25  68 

17:55:12.7 

69 

1.4 

32 

207 

08 

29  68 

19:51:30.1 

72 

1.4 

34 

262 

09 

01  68 

05:59:33.3 

75 

0.8 

19 

160 

08 

15  69 

17:45:52.7 

23 

1.9 

25 

24 

NO  CONVERGENCE 


TABLE  X 


Rei 

ocation 

of  Events  with 

Origin 

Times  Restrained 

Mo 

Date 

Shift 

Day  Yr 

Origin  Time 

Depth 

SDV(R) 

KM 

AZ 

10 

06 

64 

N/A 

09 

11 

66 

N/A 

09 

26 

66 

N/A 

07 

07 

67 

N/A 

06 

28 

68 

20:35:01 

82 

1.4 

21 

249 

06 

30 

68 

05:04:14 

55 

0.7 

7 

291 

07 

01 

68 

03:11:19 

126 

2.2 

42 

186 

07 

04 

68 

06:46:04 

68 

1.1 

18 

306 

07 

13 

68 

06:06:00 

79 

1.3 

24 

209 

07 

14 

68 

18:12:48 

66 

1.9 

21 

267 

07 

16 

68 

22:23:13 

83 

1.5 

33 

207 

07 

19 

68 

18:49:08 

84 

1.8 

43 

239 

07 

23 

68 

20:51:55 

76 

1.3 

36 

239 

07 

25 

68 

N/A 

07 

26 

68 

12:44:07 

20 

4.3 

54 

323 

08 

23 

68 

12:01:25 

78 

3.5 

39 

225 

08 

25 

68 

17:55:08 

24 

0.8 

8 

203 

08 

29 

68 

19:51:24 

16 

1.1 

29 

316 

09 

01 

68 

N/A 

09 

03 

68 

17:46:00 

91 

2.0 

20 

218 

08 

15 

69 

07:15:46 

56 

3.5 

76 

293 
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TABLE  XI 

Relocation  of  Events  with  Origin  Times 
Restrained  and  Travel  Time  Corrections 


RELATIVE  LOCATION  STUDIES 


Although  the  exact  epicentral  location  of  the  events  studies  do  not  have, 
m general,  any  bearing  on  the  discrimination  problem,  knowledge  of  these 
locations  could  be  helpful  in  relating  the  events  to  local  tectonic  patterns. 
Of  special  interest  are  the  clustered  events  situated  around  30°N  and  95°E. 
These  events  are  very  similar  in  their  characteristics  and  occur  in  a limited 
area.  The  epicentral  locations  given  by  NOS  are  given  to  an  accuracy  of  .1° 
and  are  scattered  over  an  area  of  50  x 50  km.  Using  a different  network  for 
each  event  would  result  in  a different  network  bias  in  location  of  each  event. 
Thus  even  if  all  events  occurred  at  the  same  point,  their  computed  location 
would  scatter  over  a wide  area.  By  using  common  networks  for  several  events 
it  can  be  expected  that  their  relative  locations  can  be  determined  with  higher 
accuracy,  although  the  absolute  locations  could  be  wrong.  Therefore  the 
following  approach  was  adopted.  Common  networks  were  chosen  for  pairs  of 
events  and  the  differences  in  latitude  and  longitude  (AA±J ) were  computed 
for  each  pair.  The  values  obtained  naturally  do  not  cor-'fitute  an  internally 
consistent  set;  therefore,  the  values  are  adjusted  using  least  squares. 
However,  this  procedure  is  not  quite  correct,  since  the  errors  of  the  various 
AAij  and  A<t>ij  are  not  independent,  because  they  were  computed  by  using  over- 
lapping sets  of  observing  stations.  The  correct  procedure  to  use  is  the 
joint  epicenter  determination  method  of  Douglas  (1967);  see  also  Ahner, 
Blandford,  and  Shumway  (1971)  for  which  no  working  computer  program  is 
available  at  present  at  SDAC.  Nevertheless  the  applied  procedure  should  also 
remove  the  network  biases  and  result  in  acceptable  relative  location  provided 
that  the  travel  time  anomalies  behave  consistently  over  the  area  studies.  The 
first  atcempt  to  use  the  method  outlined  above  resulted  in  an  average  rms 
inconsistency  of  relative  locations  larger  than  .1°.  This  was  judged  to  be 
too  large,  so  combinations  with  large  travel  time  residuals  were  discarded 
and  the  procedure  repeated.  The  rms  error  decreased  to  .07°  (which  is  still 
quite  large)  and  the  locations,  although  different  from  those  given  by  NOS, 
did  not  seem  to  occupy  an  area  of  different  size.  The  relative  locations  of 
the  same  events  are  shown  in  Figure  16a,  together  with  locations  by  NOS  in 
Figure  16b.  Both  depth  and  epicenter  calculations  thus  show  that  discrepancies 
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Figure  16a.  Recomputed  locations  of  events  around  30°N-95°E 


c're  present  in  the  travel  time  readings  and  the  travel  time  anomalies  which 
prevent  more  accurate  determination  of  depths  and  epicenters.  The  only  posi- 
tive statement  that  can  be  made  is  that  the  events  do  not  occupy  a wider  area 
than  originally  supposed,  and  there  are  no  indication  that  they  are  confined 
to  a considerably  smaller  area,  a narrow  fault  trace,  say.  Since  event  2 
(in  Table  1)  had  only  3 P readings  it  was  omitted  from  the  calculations. 

Summarizing  the  recomputation  of  event  locations  and  depths,  it  can  be 
concluded  that  all  attempts  to  refine  locations  and  depths  were  defeated  by 
poor  data  quality  and  no  significant  and  demonstrable  improvement  was  achieved 
The  calculations  indicate,  however,  that  there  is  no  reason  to  believe  the 
events  studied  are  considerably  deeper  than  the  depths  given  by  NOS,  so  it 
is  unlikely  that  the  low  Mg  values  are  due  to  great  source  depths. 

Mg-m^  Values 

Body  wave  magnitudes  were  recomputed  incorporating  data  at  teleseismic 
stations.  The  additional  magnitude  values  computed  are  given  in  Table  XII. 

The  magnitudes  derived  with  station  corrections  from  Report  //296  (Der,  1973) 
were  combined  with  m^  determinations  at  the  teleseismic  stations.  We  did 
not  attempt  to  derive  station  corrections  for  the  teleseismic  stations.  The 
“b  vdues  at  the  additional  stations  agreed  well  with  the  values  of  at  the 
network  used  in  Report  0296,  which  eliminates  the  possibility  that  there  were 
relatively  high  values  due  to  network  bias  at  near  stations.  The  new 
values  differ  little  from  those  given  in  Report  0296.  The  Mg  values  remain 
unchanged,  since  the  addition  of  new  stations  did  not  result  in  more  surface 
wave  detections.  The  new  Mg  and  values  are  given  in  Table  I.  The  depth 
correction  .008h  used  in  Report  It 296  and  based  on  NOS  depth  was  removed  from 
the  Ms  values  in  this  table,  since  in  our  opinion  the  body-wave  data  are  not 
of  sufficient  quality  to  rule  out  the  possibility  that  these  events  are  deep. 

It  then  the  events  are  possibly  shallow,  then  application  of  this  correction 
begs  the  question. 

Marshall  and  Basham  (1972)  derived  their  depths  from  amplitude  measure- 
ments of  the  dispersed  surface  wave  train  at  various  dominant  periods.  This 
was  not  possible  in  our  case  since  most  of  the  surface  wave  observed  were 
short  transients  which  could  not  be  analyzed  without  digitization  and  Fourier 
analysis.  The  quality  of  the  available  photographic  records  was  poor  and  it 
was  not  practical  to  digitize  them.  Besides,  calculations  by  von  Seggern 
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TABLE  XII 


Body-Wave  Magnitudes  (ni  ) 
at  Teleseismic  Stations 


Event  Date 

Station 

64/10/06 

COL 

4.67 

BUL 

4.48 

SHI 

4.47 

66/09/11 

BUL 

4.14 

COL 

4.03 

SHI 

4.30 

66/09/26 

CTA 

4.48 

COL 

4.59 

BUL 

4.29 

67/07/07 

BUL 

4.76 

COL 

4.89 

CTA 

4.48 

SHI 

5.01 

68/06/28 

BUL 

4.62 

COL 

4.93 

CTA 

4.90 

SHI 

4.62 

TRI 

5.17 

68/06/30 

BUL 

4.68 

COL 

4.79 

68/07/01 

COL 

4.42 

TRI 

4.90 

68/07/04 

BUL 

4.84 

COL 

4.85 

NUR 

4.83 

SHI 

4.78 

68/07/13 

COL 

4.87 

SHI 

4.67 

68/07/14 

BUL 

4.94 

SHI 

4.79 

COL 

4.98 

NUR 

4.87 

TRI 

5.26 

68/07/16 

BUL 

4.83 

COL 

4.92 

NUR 

5.08 

SHI 

4.49 

Event  Date 

Station 

68/07/19 

BUL 

4.90 

COL 

4.85 

’ UR 

4.91 

SHI 

4.53 

68/07/23 

BUI 

4.96 

COL 

4.78 

NUR 

4.83 

SHI 

4.61 

68/07/25 

BUL 

4.90 

COL 

4.98 

NUR 

4.94 

SHI 

4.78 

68/07/26 

BUL 

4.96 

COL 

4.82 

NUR 

4.80 

SHI 

4.58 

TRI 

4.91 

68/08/23 

BUL 

4.66 

COL 

4.92 

NUR 

4.94 

68/08/25 

BUL 

4.88 

COL 

4.85 

CTA 

4.66 

NUR 

4.67 

68/01/29 

BUL 

4.68 

COL 

4.92 

NUR 

5.04 

SHI 

4.64 

68/09/01 

BUL 

4.50 

COL 

4.85 

SHI 

4.60 

68/09/03 

COL 

4.82 

CTA 

4.88 

NUR 

4.91 

69/08/15 

BUL 

4.91 

COL 

5.09 

SHI 

4.78 
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using  Harkrider’s  theory  of  surface  wave  generation  showed  that  the  method 
used  by  Marshall  and  Basham  to  determine  depth  is  extremely  unreliable.  Thus 
the  separation  of  earthquakes  and  explosions  achieved  by  Marshall  and  Basham 
must  be  partially  due  to  spectral  differenc  r between  earthquakes  and  explo- 
sions. Without  the  "depth  correction"  the  separation  of  explosion  and  earth- 
quake is  less  pronounced,  and  low  Mg  earthquakes  are  liable  to  be  mistaken 
for  explosions. 


This  problem  is  especially  severe  for  relatively  small  events  where  the 
surface  waves  can  be  seen  only  at  nearby  stations  at  low  signal-to-noise 
ratios.  This  is  precisely  the  case  we  are  dealing  with.  Without  the  depth 
correction  most  of  the  events  studied  here  seem  to  have  an  anomalously  low 
Mg  compared  to  n^,  and  could  potentially  be  mistaken  for  explosions  if  the 
Ms_mb  discriminant  is  used. 


.Effect  of  Abnormally  Thick  Crust  on  Surface  Wave* 

The  source  region  studied  is  anomalous  in  its  crustal  thickness.  Estimates 
of  crustal  thickness  from  surface  wave  studies  are  in  the  range  70-90  km 
(Gupta  and  Narain  1967).  The  anomalously  thick  crust  presumably  results  from 
the  under thrusting  of  tne  Indian  subcontinent  under  Eurasia.  To  ascertain 
that  the  low  surface  wave  excitation  is  not  due  to  abnormally  thick  crust  we 
have  calculated  the  surface  wave  excitation  factor  ^ (Harkrider,  1964). 

From  Harkrider  (1964)  and  from  equations  (12)  and  (24)  in  Ben-Menahem 
and  Harkrider  (1964)  one  may  derive  the  equation 

„ 1/2  w 

R o 

where  M(u)  is  the  source  strength  (M(o)  = moment),  CR  is  the  Rayleigh  phase 
velocity,  and  Wg/WQ  is  the  ratio  of  model  amplitudes  at  source  depth.  Note 

that  Wg/WQ  is  unity  for  surface  sources,  and  greater  at  any  depth  in  a thick 

crust  than  in  a shallow  one  because  the  fundamental  mode  anplitude  decays  more 
slowly  if  the  velocity  increases  slowly  with  depth.  CR  decreases  as  the  crust 
increases^  thickness  and  von  Seggern  (1971)  has  shown  that  ^ is  proportional 
to  (l’CR)  (where  U is  the  group  velocity)  which  also  would  be  expected  to 

increase  in  a thick  crust.  Thus  one  would  expect  that  the  effect  of  a 70-90  km 
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crust  would  be  to  Increase  Mg  if  measured  on  the  thick  crust,  not  to  decrease 

it.  And  in  fact  explicit  calculations  give  increases  of  .2  to  .4  M units 

s 

over  a Gutenberg  crust. 

Now  let.  us  suppose  that  the  amplitude  measurement  takes  place  on  thin 
crust  and  that  the  wave  has  crossed  via  a smooth  transition.  It  seems  plausible 
that  the  amplitude  would  revere  to  what  it  would  nave  been  for  a source  in 
the  thin  crust.  In  fact,  McGarr  (1969)  obtains  theoretically  and  experimentally 
a shift  of  magnitude  not  significantly  different  from  U~3^2  for  an  ocean 
continent  transition.  Thus  for  measurement  on  a thin  crust,  we  would  expect 
t' e net  effect  to  approximately  vanish.  Note  that  this  failure  to  explain 
low  surface  waves  by  a completely  elastic  mechanism  does  not  rule  out  the  low 
Q proposed  by  Lambert  (1974)  as  an  explanation. 

The  dominant  period  of  Rayleigh  waves  observed  is  10-15  sec  which  indi- 
cates that  most  events  considered  are  shallow  since  greater  source  depth  would 
drastically  decrease  the  amplitude  of  these  waves  relative  to  the  longer 
periods  which  we  do  not  observe. 


SEARCH  FOR  ADDITIONAL  ANOMALOUS  EVENTS 


Search  for  low  (Mg-n^)  events  was  performed  in  a wider  area  around  the 
region  previously  studied.  All  NOS  events  occurring  between  25-45°N  and  68-100°E 
in  the  time  interval  June  4,  1968,  and  September  30,  1972  were  used,  except  for 
those  already  discussed  in  SDL  296,  and  those  which  had  no  listed  by  NOS. 

The  time  interval  given  was  determined  by  the  availability  of  the  data  from  the 
WWSSN  station  Kabul  (KBL)  , Afganistan.  This  station  operates  at  very  high 
magnifications  on  both  the  short-period  and  long-period  instruments.  It  is 
located  inside  the  region  searched,  and  is  therefore  quite  suitable  for  detecting 
the  events  to  be  studied.  It  can  be  expected  that  the  single-station  values 
would  be  unreliable  since  path  and  station  effects  for  short-period  P waves  can 
result  in  standard  deviations  of  the  order  of  .5  magnitude  units.  Single-station 
Ms  values  have  a smaller  scatter  due  to  path  and  station  effects  but  can  be 
influenced  by  radiation  patterns.  In  view  of  these  factors,  the  following 
strategy  was  adopted:  values  as  listed  by  NOS  were  used  (which  were 

presumably  computed  by  averaging  values  of  several  stations)  jointly  with 
Ms  values  determined  at  KBL  alone.  Previous  experience,  Der  (1973),  showed 
that  values  given  by  NOS  are  in  general  higher  than  these  determined  by 
SDAC.  Too-high  values  would  minimize  the  chance  of  missing  some  truly 
anomalous  events.  Similarly  some  Mg  values  will  be  too  low  at  KBL  due  to 
radiation  patterns.  These  events,  if  truly  normal  in  their  charac- 

teristics, can  be  eliminated  later  after  closer  examination.  Thus  the 
combination  of  NOS  and  KBL  stations  Mg  values  can  be  used  to  sift  through 
the  initial  data  set  and  eliminate  most  events  with  normal  Mg-n^  values. 

Table  XIII  gives  the  events  in  the  initial  data  set  with  the  Mg  and  values 
used.  Events  which  are  not  given  here  but  listed  by  NOS  were  eliminated 
because  of  undetectable  signal  amplitudes,  some  operational  difficulty  at 
KBL,  or  interference  with  other  events.  It  is  of  course  possible  that  some  of 
the  events  eliminated  because  no  LR  ’/as  detectable  were  truly  the  events  we 
were  looking  for.  The  Mg  values  were  computed  by  Marshall  and  Basham's 
method  using  depth  corrections  with  depth  taken  from  NOS.  Figure  17  shows 
the  plot  of  Mg-n^  values  obtained.  In  all,  165  Mg-n^  pairs  are  shown.  The 
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TABLE  XIII  (Continued) 

Body  and  Surface  Wave  Magnitudes  at 
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TABLE  XIII  (Continued) 

Body  and  Surface  Wave  Magnitudes  at 
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Body  and  Surface  Wave  Magnitudes  at 
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no  Ms  - 1.38m,,  - 3.45  was  used  as  the  dividing  li„e  between  "anomalous"  and 
normal  events;  it  lies  midway  between  the  earthquake  and  explosion  lines  of 
Marshall  and  Basham  (1972).  Eleven  points  lie  below  the  line  (the  events  in 
question  are  marked  in  Table  XIII  by  a star).  One  of  them,  the  event  in 
Sinkiang,  is  a presumed  explosion. 

To  investigate  these  eleven  events  further,  film  chips  were  ordered  from 
NOS  for  the  stations  NDR,  TRI,  SHI,  POO,  SUL,  QUE,  BUL,  ME,  CTA,  COL,  NDI  in 
addition  to  KBL.  P wave  and  Rayleigh  wave  amplitudes  were  read  and  magnitudes 
computed.  The  recomputed  magnitudes  are  shown  in  Table  XIV  and  are  plotted 
against  Marshall  and  Basham's  95%  confidence  limits  for  earthquakes  and 
explosions  in  Figure  18.  The  filled  circles  in  this  figure  are  the  M -m 
pairs.  The  figure  shows  that  only  two  of  the  points  lie  close  to  the°95* 
confidence  limit  of  Eurasian  explosions;  one  of  these  points  is  an  explosion. 
The  remaining  points  are  close  to  the  95%  line  for  earthquakes.  Since  there 
is  a considerable  uncertainity  in  determining  these  lines,  it  is  not  clear 
that  anything  is  unusual,  since  the  number  of  "anomalous"  points  is  about  5% 
of  tae  total  events  considered.  Thus  the  population  searched  is  not  much 
different  from  that  presented  by  Marshall  and  Basham.  The  most  anomalous 
earthquake  is  located  in  the  Tadzhik  SSSR,  but  the  point  is  based  on  only 
three  m^  and  one  reading. 
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SUMMARY  AND  CONCLUSIONS 

The  analysis  of  low  M ovents  In  w-tm  i 
„ , s tS  in  the  Himalayan  region  using  an  expanded 

network  gave  the  following  conclusions: 

LJhe  Val"S  "*  COnflrmCd  “ be  "latlVely  l0“  C0"pared  “ \ values. 
Additional  P-wave  mleselsmlc  leadings  did  not  change  the  ^ values  apprecianly. 

2.  Clear  dilatations  were  observed  at  for  many  events,  indicating  that 
they  are  earthquakes  rather  than  explosions. 

<•  Depth  and  location  calculations  using  various  constraints  and  e common 
network  failed  to  provide  refinement  in  the  source  location  or  depth  due 
to  the  poor  quality  and  internal  inconsistencies  of  data.  The  calculations 
indicate,  however,  that  most  of  the  events  are  shallower  than  80  km. 

4.  PP  readings  at  teleseismlc  stations,  although  doubtful,  also  indicate 
that  most  of  the  events  are  shallow  or  normal  1„  depth. 

5.  The  long  period  S to  Rayleigh  amplitude  ratios  classify  the  events  as 
earthquakes . 

6.  Anomalously  thick  crust  in  the  source  region  is  not  likely  to  cause 

the  low  M values, 
s 

7.  The  event  cluster  centered  at  30"N  latitude  and  95°E  longitude  consists 

of  many  similar  events  with  almost  identical  long  period  seismograms,  hove 
waves,  If  present,  are  unclear. 

8.  Three  events  on  the  Himalayan  front  show  long-period  S,  two  long-peri„d 

and  two  hove  waves.  The  long  period  S to  Rayleigh  amplitude  ratio  is 
characteristic  of  earthquakes. 

9.  The  short-period  S to  P ratio,  where  available,  is  characteristic  of 
earthquakes . 

10.  Search  for  additional  anomalous  events  in  a wider  area  did  not  indicate 
any  additional  contiguous  areas  where  such  events  occur.  The  low  M events 

found  can  be  explained  by  normal  scatter  of  data  and  amount  to  a f«  percent 
of  the  total  data  studied. 
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The  above  results  suggest  that  even  low  (l^-n^)  events  such  as  those 
discussed  by  Der  (1973)  can  be  shown  to  be  earthquakes  by  careful  analysis 
of  the  short-period  data  from  close-in  stations.  Long-period  shear  to  Rayleigh 
amplitude  ratios  seem  also  to  classify  them  as  earthquakes.  The  data  quality 
is  poor  enough,  however,  that  a fully  satisfying  discussion  of  the  events  is 
impossible.  In  particular  we  remain  uncertain  about  the  cause  of  the  low 
(Ms‘mb)  values;  and  we  cannot  be  certain  that  the  long-period  shear  wave 
discriminants  are  completely  reliable  because  the  measurements  were  made  at 
distances  oc  only  5°  for  which  distance  there  are  no  data  in  the  discrimination 
data  bases  of  von  Seggern  (1972)  or  of  Blandford  and  Clark  (1974).  Even  if 
such  data  were  available  it  would  have  to  be  carefully  analyzed  from  a 
theoretical  point  of  view  to  ensure  that  the  unusually  thick  crust  in  Tibet 
does  not  bias  the  results.  These  remarks  point  out  possible  subjects  for 
future  research,  which  should  become  possible  with  the  aid  of  high-quality 
data  from  the  expanded  VELA  network. 
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